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Processing Seismic Data 

■ c data and in particular, to processing 
tte pres e„t invention rela.es ,o P— -™ ^ or iemove te swell 

nrarine seismic data that rnc.udes swell no.se so as 



noise. 



■. emitted from a source and is detected by a sensor 
b . seism, anrvey, acoustic energy . . * ^ fc ^ Some 

(generally referred ,o as a "recerver , ,oc«d ^ ^ ^ ^ by 

-* — — 1~ ^rln abou, me g eolo,oal stincbare of tire 
^^^edtm tire reflected acoustic energy incdeo, on tire 



earth's interior can 
receiver. 



- ™<\ receiver are disposed within a 
b . manne aeismrc survey bora — ^ ^ —n of „„e marine 
water column, such as the sea. ng» js OTSpended in 

seismic survey in which an aeon* energy source^ ^ ^ f ^ ^ „ 
fc water column 2 from a survey vesae ^ ^ ^ ^ , „ 

.ample, hydrophonea, * * » *^ , Such . reeerver 

array is known generally aa a rowe 

, is aomated >. emhs acoustic energy. Acoustic energy 
When tire aoouatic source 1 - ^ fce water colulM , 

emitled by .he source prepay ^ ^ dQWnwards 

passes into the earth's rnteno and ^ ^ ^ 

throug h tire earth's interior until it nndergoea par. ^ 

• , • . Ppflected acoustic energy p^p s 
wim in .he earth's m.ertor. Reflected ^ ^ ^ ^ M , he 

fcorrgh me earth's interior 1 and throng of , he 

— 5 s rir nrs « — «■ - * - 

acoustic source 1 are — ^ be ^ in . processor S located on .he 
eUcttical or optical h* Th * ^ ^ principl e of aeramre 

second survey vessel or may be stor 



m&c Folio No P52H4GB nrocessing of the 

■ ^„t the earth' s interior by suitable processing 
surveying is to obtain information about the earth 



seismic data acquired at the receivers 5. 



, Wfl the ear th's interior acts as a partial 

2£- - • — 5 — ■ 

refl ec.or for acoustic energy. The rf ^ ^ 

nation of tite acoustic source ^ ^ ^ ^ M 

geological feature within the earth tna. 

. martae seismic data is that data acquired in rough sea 
One proh,em » processiug towed mann. ^ „ ^ « 4 

Editions may contain swell no.se. « ^ gmeiaKd by . cough sea 

occurs at low frequencies, generally * „ Kismic data (i.e., the events 

in ,ecferes destine,, «* - *"* ^ features wtt dn the earth). Swell 

fxon, reflechon of setsnnc energy ^ 

.ountherofpnorartntethodsforanenuationofsweunorseare.nown. 

to "An.onta.ic editing of noisy seismic data", 
OA. McMechan has proposed >n ^ ^ fa ^ ^ ,„„ 

Geophysics! Prospering, Vol. 37, « ^ ^ ^ ttaC e is discarded if its low 

frequency energy content of a se,snu jiu S method has the severe 

creeds a given tnrestioiu. 
frequency energy content exceed, g ^ 
diaadvartage ft., tire desired stgnaltnfltetiacets 

d Robert s Wo,slaw have proposed, in "Spatially weighted 
„.Truntan Holcomhe and Rohert S. J ^ ^ of Expl0ratl0n 

* ^ 2 rr C. «) n .!> involves locating awe,, noise in a 
Oeophyaieiats, ppl 157-1 160 (« ^ mtHpolatio „ whic h makes 

tra oe and replacing ,he affected pontonof *e ^ ^ ^ ^ 
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♦„ «rr,re<5s Furthermore the aip =>g 
costly and often inaccurate process. 



dominant dip. 



a in "Can we surgically remove 

swdl n „ iS e -t: ; ^P,, doma, .» - - 

to, is designed and apphed m » x ( * ^ noise . nol 

efficl en, a, filtering ou, swefi no,se when ^„ ^ 

ellffli „ aK d by * — cann. — — - - 

fitter* a coherem enhancement process an 



from a coherent signal 



... ...dative Seismic Noise and 

AU o** has proposed, in a paten, ^ fer attenuadng 

Terence Attenuation Mednod" ^ ^ M m e,hod 

ta various types of coherent no.se — e ^ ^ 

applies „ adaptive (,-x) *» space ffl er m P ^ ^ ^ ^ 

that contain me «•"■ not statl0 n W are adapdvely 

Terence components ma, are ^ J ^ 
fitted Thrsmemod can be qurteeffecfivebu 

, a method of attenuating swell noise, in One SGE if 
LC-K »» *> proP ° Sed 3 ™* Thrs is a "one step ahead" wiener predion 
^ ttttemationa, Meeting ^ .„ Co „ ermt ^ and 

filtering proofs earned ont m the U (fteq 



Cerent swell noise are bom enhanced. 



• ninFU«eringinmeF-XDomain».SEG6 3 -An„na. 
S.Sp te has proposed, in "Adapdve Dtp FtKeng ^ t0 

agonal Meeting (« • "^^P " " " 

a „ „ dips widdn me da,a are — ^^ tade , w «n the hrpu, data are 
pK dicfion filter. All the dipping «* ^ ffi all *e dips need 

L solved for using a, east square approaches, 



to be accounted for. 
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inv0 We comparing • - ^ w thre shold vahue. Xo.se above 

for to presence of noise, usuaUy „ ywy smsWv e • «. 

*. specified to- * «-»-■ S " Ch J ^ te 6ala deluding some of .he 

te esho,d, and - — »» « ^ 
sigMl of tate^s, or ,o leave some res.doalno.e. 



signal oi uitvx^-- 

• a*** the data being in 
, j ^rnre^sine seismic data, we u<u 

*. freuuency-space *— • «" fot fte ,„ OTCy a, which swell no.se ,s no. 

prei icfion fiber from fire se,sm,c ^ ( at a sec o„d fie,uency a. whach 

swell noise is presem, j 
frequency. 

• <r nm the signal of interest in the data, and 
^ fiction filter separates the swell noise £££ ^ ^ 

---rr^ - swei1 noise without 

provides effective attenuation of both 
affectingthe signal of interest. 

m . . method of marine seismic surveying 
A SK ond aspec, of fire — acquiring sersmfc dala ,n ,e 

^ — ■ «^ " da , a accordtag ,„ a nae,od as defined above, 
frequency domain; and processing 

^tns for processing seismic data, the 
A1 « aspec. of - — p.ovidpjP^J tos ^ _ foI 

dala being in -he ^ at a « frequency a, which swell 

det ermin»,g a predion finer »- fflter 10 selsmia da,a a. , 

nnlse is no. presen.; and — •- > ^ ^ ^ norse in fiae 

.econdfiequency a. which swell no.se ,s present 



seismic da.a a. the second frequency. 
Tfce apparatus may contain a programmable 



data processor. 



M&C Folio No P52144GB crtlirrP of 

A n „ arrangement comprising: a source 01 
. rtfcP invention seismic surveying arrangeme r 
A fourth aspect of the invention freq uency-space domain; and 

„e fnr ar auiring seismic data in me 
acoustic energy; means for acquiring 

an apparatus as defined above. 

• mvides a storage medium comprising a program for a 
A fifth aspect of the invention provides rf ^ ^ 

data processor of an apparatus as define n aho. J ^ ^ ^ 

provides a storage medium containing a program 
perform a method as defined above. 

, • ♦:„„ are defined in the remaining claims. 
Preferred features of the invention are define 

• ti«n will now be described in detail, by way of an 
A preferred embodiment of the invention wiU now ^ 
illustt ative example, with reference to the accompanying figur 

• ^ fr^nuencv domain, obtained by the 
. Kmiral seismic data trace, in the frequency u 
Figure 2 illustrates a typical seisms 

seismic survey of figure 1; 
invention; 

seismic data; 

W fc block diagram of an apparatus according to the present 
Figure 6 is a schematic block cnagran 

invention. 

• A * one of the receivers 5 in the seismic surveying 

:::::: i^^-.-*— 
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prM s« o, a componen, -« ■ "J ;° Ke „ aces ipCude swell noise, 

eismio data .races are shown » F, g urc^ - 5 ^ 4(a) md 5(a , aft. 

^ „ Figures 4(h) - *0» «** ~ * „ ttme , . 0, so the. a trace 

removal of swell noise). The acousttc omce . ^ — q{ fte 

„ -he variation in **"*«*~> - *> ** *° * "t 
aM us,ic source. » will he seen ftom ^ „„ slgna , a, fire 

^ „ acmarod, mere is in."* .. * ^ ftom ^ energy 

re eei,er. The firs, even, rn .he scsmtc d.te $ ^ ^ ^ ^ 

« has ftaveiico dtrect ftom Ihe — -J * ^ ^ actuaD on of the 

* "Tl ,? to o» as .he "direct evenO depends on the 

acoustic source and this even, (winch .s to ^ ^ ^ distance 

sp eed of propagation of acoustic energy through ^ ^ ^ by 

the acoustic source and tire waver. The dn ^ 
be tween me ^ ^ enterei to ear* a urten 

„ receiver will generally sample the 
however, a receiver e intervals of 

„ f, (-ha. ., *e paranrcter „ a»pled a ^ ^ ^ ^ rf , ^ 

,„«„». « nray he, for ^ s „ pUn6 opeta ,i„n. 
of discrete values for tire parameter, on 

• „v ha Fourier transformed to produce a 

n. seismic date Pace in the time dom^m ^ ^ Rgure 2 „ 

record of me amplitude of me »- P* * domata , as would he 

. schema 01— of a typical date ha ^ ^ ^ ^ „. 
ohteined hy transforming one of .he date ^ ^ Mce . tha 

spaM domein hy, for ex»P>e, « ^ ^ t ^ 

L douteh, is ohteined hy a saanphng ^~ „ of f w where 

formed date trace of figure * •— «*•' torem . The date trace is 

f „„ , s re ,a,ed to .he sampfing fi^ * ^ ^ ^ . 

shown as a continuous nac. m F.gure f § ^ of discrete valu es 



K! "" 3B . ,„ ftMuenc v-space domain of figure 2 

,„ ,han a continuous trace. The data trace m the freq»«n P dex 
rather than a com, ^ 6equencies tj , where 

W» *"*" *° ^ ,° , ; wteh fj does no. exceed the Nyouist frequency, 
j extends to the highest value J at which 

• ,.„ havimt a maximum fiequency of around 
SweU noise is a, ow^nenoy noise,, yp,« f ^ wo 

30 Hz. The ^^T^ -fit «te srgna, of mterea, - -* 
reg i„„s - a low fieoueney re^on 9 « ^ rf toterest . „ 

no- and a higher — ' ^ ^ nols e and the of — 

fi6urc i the « IT1 « * - - * " - " 

tv,P data at frequencies below il 
iff. is not zero, the oam ai i 



If f L is not zero 
noise other than swell noise. 



• al orediotion filter is detennined fiom me seismio 
„ .he memod of the — . a stgna ^» ^ ^ _ „ js> at . 

- - — * r r ri * - a PP ,ied . d> : « 

frequency of fm, or greater. * P ^frequency between t and 

_d frequency a, which "^"^ ^ „ secon d frequency. * 
fH . s „ as to attenuate swel no, « ^ fcr ^ „ more ote 

melh od is preferably repeated ,o ^ ^ for example by deriving 

fre q „encies a, which awe,, no.se is P-j " ^ ^ ^ md ^ .his 

. fine, fiom dam a, a third ^«y(^h ^ ^ i ^ ^ ( „ 

fiction filter to ^ ^ prrfOTed ^dtment, an assocated 

which swel, noise ta present - ^ ^ ^ ^ ^ ^ bom k * f „, 

range. 

„, ,he idea, ease, the prediction fi„er .hot was « 
den ved from seismic dam for that frequency. In P 
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crit erion would ^ « <■» 1 from the seismic da« for ft* 

frequency in the <*» fL » » ™° U " Mon fflKt from seiamic data .ha, 

„ - hu, — — 7' ^ 0 t res , According .o *• — • 
cmtate swell - * - ^ :; Mt , > den ved from me scsmic 

before, ft. P red.c.,on Mr for a » J fc 6equency rmge wl ftm 
dat a * a fluency «ha, U -« " * ^ app— n, Prided tha, the 

wWc h swe„ noise «-J I - ^ ^ similaI , tor the frequency < 

d.aracterirtcaoffhes.gaalofuiterestar 



aim xv* - — 

, t f r a single frequency at 

„ principle, a single P— « * "** ^ 

wbi O swell noise is no, present could * * ^ „ „, fiUeI is 

nois e is presen, This howeve, has ^ ^ ^ be „ oUy at 
for any reason, then - * ^ ftat . plummy of 

ever y ^ency ,n dre swell noise ^ „ e used so ta , 

pr *,,n fhers, each fi „ ^ prefro h,e .a, £ 

freq ue„c ie s at whi ch s ^ ^ ^ ^ 

preferably determined at a frequen y 

♦t, filter to be applied to the 
b . p« Prefer em— of * ~ * ^ „ „ is, as 

da ,a a, fluency *. «■ d — 7 „ ise occ «, and is ft. .owes, 

not ed ahove, ou,side ,he range a, winch s- n ^ - ^ 

outsi de fte range a, whioh swell no,se occu*. Th gn 
... • . -I- thnse for frequency tj,- 



outside tne rmi B ^ *» 

fol should he similar ,od,ose for frequency f. 



, f is derived from the seismic 
ta th , s emhodiment, the fdter to - — - - 

a, fa. Again, the ^ ' S 
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noise occurs, and the signal of interest at frequency should be similar to the signal 
of interest at fL+i- 

Furthermore, in this embodiment, for each frequency within the frequency band for 
which swell noise occurs, the associated prediction filter is determined from the seismic 
data at a frequency fi-L+H+i- Thus, finally, the prediction filter for the highest frequency 
fw in the swell noise band is derived from the seismic data at frequency f2H-L+i- 

The frequency range in which the prediction filter is derived thus starts at fn+i, and this 
frequency fulfils two important criteria. Firstly, data acquired at this frequency, or at a 
higher frequency, is assumed to be free of swell noise. Secondly, this frequency is as 
close as possible to frequency ft, so that the signal of interest should be similar at 
frequencies ft+i and ft. Furthermore, in this embodiment each filter is applied only to 
data at one frequency, so that the effect of a corrupted filter is minimised. 

Figure 3 is a block flow diagram of one embodiment of the present invention. 

Initially at step 1 seismic data are obtained, for example, using a towed marine seismic 
surveying arrangement of the type shown in figure 1 . The seismic data will consist of a 
number of traces generally similar to those shown in figure 4(a) or 5(a). 

The invention may alternatively be applied to pre-existing seismic data. In this case, 
step 1 is replaced by step 2 of retrieving pre-existing seismic data from storage. 

At step 3, the seismic data are Fourier transformed to produce seismic data traces of the 
general type shown in figure 2 that contain information as to the frequency dependence 
of the amplitude of a seismic parameter. 

It should be noted that, in an embodiment where stored seismic data are retrieved from 
storage in step 2, step 3 may be omitted if the retrieved seismic data are already in the 
frequency domain; this is indicated by a broken line. 



10 

M&C Folio No P52 144GB 



M&C Folio No P52 144GB .-.j f or 

i ♦ a filter F, is detenmnea ior 
• ; e initialised and at step t> a nner r s 
At sBp 4, a frequency count, . ^ ^ ^ ^ ^ h (he seismic 

„„ initial frequency compouent f, m the r g ^ at 

d ». Conveniently, step 5 , initially performed or he * , ^ ^ 

* *v»ot i« for frequency U - ±l 
which Wl . noise rs present tha * ^ ^ ^ embo6ime „ t 

nave set ft. counter i equal to I}. ^ ^ ft ^ fm . 

,he prediction filter for fieqnency ft rs detemun 

f for all traces, are operated on by the filter Fs 
At step 6, the data at frequency C for all t« , 

«. fin step 5 according to equation 1. 

determined for fiequency f, in step o, 

• a whether the filter has b«n applied a pre-determined number 
At step 7, it is determined wh .her to ^ osH for processl „ g a 

„ times. The predetermined nu« J ^ g ^ ^ R „ „. 
pa ,icular dafa so, If step 7 yrelds no de ^ ^ ^ ^ , f ^ „ yK „ 
nppiied ,o to resnba ** rf m ^ 8 ^ 7 1 repeated nnfil a determination 
netermination is still not obtained, steps 8 and ar 



is obtained. 



• +rt w nnnlied three times (n - 3), in 
For ^ a .or could specify - — - to the 
which case step 8 womd be canaed u.^ ^ ^ ^ method 

* - -r z ~ - - — - — of ,,mes ,he H,er 

ma y compnse the step <nui a 



should be applied. 



k . men, s.ep 7 is replaced by a ,ual«y oonhol sfep (no. show.), 
m an al.ema.ive embodunent, s ep 7 P ^ ^ ^ ^ , pfe . 

M t aU t:":"e residua, ob,,ned an one sfep may be 
determined number of times. 

compared with the previous residual. 

■ „ whether die counter indicates to. the ptoceas has been 
A, step 9, is determtned wheto. ■< ^ ^ ^ lf , hevalueot 

cam ed out for a„ frequences . 6 » ^ ^ ? , „ no „ deBnninati0 „, 

ft. counter i is no. equal to H so f, s ^ ^ ^ fcr , he 

,he frequency counter i is incremented a. step 10, 
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next frequency - in this example, f L +i. Assuming that a "no" determination is still 
obtained at step 9, the counter is again incremented at step 10, and steps 5-9 are 
repeated until a "yes" determination is obtained in step 9. A "yes" determination will be 
obtained when a method has been carried out in all frequencies in the range fL to - 
that is, for all frequencies at which the data contains swell noise. , 

4 

I 
1 

Alternatively, the method may comprises performing steps 5 to 8 initially for f H , in 
which case step 10 comprises decrementing the counter. In this variation a "yes" 
determination is obtained at step 9 once the frequency has been reduced to f L . 

Once a "yes" determination is obtained at step 9, the results of the method may be 
output if desired (not shown). The results output may include, for each original data 
trace, a corrected version of the data trace from which swell noise has been removed. 
The corrected data trace may be obtained, for example, by subtracting the final residual 
obtained at the end of the iteration step process 7 and 8 from the original data trace. 
Additionally or alternatively, the results output may include the residual obtained from 
each data trace. 

Additionally or alternatively the corrected data traces may be subjected to further 
processing steps, to obtain information about the geological structure of the earth's 
interior, and this is indicated at step 12. If necessary, the corrected data traces are 
transformed back to the time domain at step 1 1 before step 12 is applied. These further 
data processing steps are conventional, and will not be described here further. 

Figures 4(a) to 4(c) and 5(a) to 5(c) illustrate results of the present invention. Figures 
4(a) to 4(c) illustrate the invention as applied to synthetic seismic data, and figures 5(a) 
to 5(c) illustrate the invention as applied to actual marine seismic data. 

Figure 4(a) shows synthetic marine seismic data obtained using a simulation method 
that generates swell noise. The data were simulated for a seismic surveying 
arrangement in which the direct arrival occurs approximately 0.8 seconds after the 
acoustic source has been actuated. Figure 4(a) shows approximately 75 simulated data 
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traces, in the time domain. Each trace shows the amplitude (horizontal axis) of the 
pressure as a function of the time (vertical axis) after the actuation of the acoustic 
source. It will be noted that the data traces contain significant variations in amplitude at 
times before the direct arrival. 

Figure 4(b) shows the traces of figure 4(a) after swell noise has been attenuated using a 
method of the invention, and figure 4(c) shows the residual which represents swell 
noise. (It should be noted that although the method of the invention is carried out in the 
frequency domain, the results for the corrected data traces and the swell noise have been 
transformed back to the time domain to enable easy comparison with the original data 
traces.) 

■ 

If figure 4(b) is examined, it will be seen that there is very little structure in the seismic 
data traces before the direct arrival occurs at approximately 0.8 seconds after the 
actuation of the acoustic source. This indicates that the invention is effective in 
removing swell noise since, in the absence of noise, the receivers should not detect any 
acoustic energy from the source before the direct arrival. It will also be noted that there 
is significantly less structure in the seismic data traces for times greater than 
approximately 3 seconds after the actuation of the acoustic source. This indicates that 
most of the structure seen in the original data traces of figure 4(a) at 3 seconds or 
greater after the actuation of the acoustic source does not result from reflection by a 

■ 

geological structure but arises from swell noise. 

It will further be noted that the events in the seismic data traces occurring at times from 
approximately 1.6 to 3 seconds after the actuation of the seismic source are much 
clearer in the corrected seismic data traces of figure 4(b) than in the original data traces 
of figure 4(a). This indicates that the attenuation of swell noise has been effective 
throughout the entire data trace. As a result, information about the structure of the 
earth's interior derived from analysis of these events should be more accurate than if the 
original uncorrected data traces of figure 4(a) had been used. 
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Figures 5(a) to 5(c) again show results of the present invention, but applied to actual 
seismic data. These seismic data were acquired using a marine seismic surveying 
arrangement in which the direct arrival occurs approximately 1.8 seconds after the 
actuation of the acoustic source. The figures show data for approximately 70 receivers, 
and these are arranged as a linear array extending generally away from the source. As a 
result, the arrival time of the direct event increases gradually from one trace to the next, 
as a consequence of the increasing distance between the source and the receiver used to 
acquire a particular trace. 

As is clearly seen in figure 5(a), the acquired seismic data contain significant amplitude 
variations for times before the arrival of the direct event and, as explained above, these 
amplitude variations arise primarily from swell noise. 

Figure 5(b) shows the traces of figures 5(a) after attenuation of swell noise using a 
method of the present invention. It will be seen that the amplitude variations in the 
traces for times before the arrival of the direct event have been significantly reduced. It 
will also be seen that many of the events in the traces for long arrival times have also 
been eliminated, and this indicates that these events arise from swell noise rather from 
reflection by a geological feature within the earth. The events of interest in the data 
traces, which generally occur at times between 1.8 and approximately 5-6 seconds after 
the actuation of the source are much clearer in figure 5(b) than in the original traces of 
5(a). 

A preferred embodiment of the filtering method will now be described. 

Let D i denote the input data. This is assumed to consist of N traces identified by the 

index i where 1 < / < N. These traces are acquired at the same time but by different 
receivers in a receiver array. Each trace will contain a signal of interest, and will also 
contain swell noise. The signal of interest in the z ,th trace will be denoted by X si and the 

swell noise value in the i th trace will be denoted by X Ni . Thus: 
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D, =X S1 +X 



Nl 



D 2 - X S2 + X N2 



D 3 =X S3 +X 



N3 



D N -X^ +X NN 

The signal prediction filter is denoted by Fs = [F Si F S2 F S3 ... F SL ], where Fs* 

denotes the k th filter component. Fs is a general term for a prediction filter derived at a 
particular frequency outside the swell noise range. The index k satisfies 0 < k < P, 
where P is the number of signal related dips. The filter length is preferably between 3 
and 10 traces. The prediction filter may be derived by any suitable technique such as, 
for example, Wiener filtering. As noted above, in a preferred embodiment the 
prediction filter for frequency ft is determined from data at frequency f H +i, the 
prediction filter for frequency fi+\ is determined from data at frequency fH+2, and so on. 

An appropriate filter Fs is applied to each seismic data trace at each frequency in the 

■ 

range (fL - fu) for all traces. The resultant signal is predictable as a function of space 
(Canales 1984, above). 

Filtering the input data with F s only will result in the following system of linear 
equations: 
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If the lower frequency limit of the swell noise range, Fl, is zero, the output will consist 
only of (a) and (b) above. 

Figure 6 is a schematic block diagram of an apparatus 11 that is able to perform a 
method according to the present invention. The apparatus 1 1 may be used in a marine 
seismic surveying arrangement in place of the conventional processor 8 of Figure 1, or 
it may be used to process pre-existing data. (If the apparatus 11 is used in a 
conventional marine seismic surveying arrangement of the type generally shown in 
Figure 1, the seismic surveying arrangement requires means for acquiring seismic data 
in the frequency-space domain and these are conveniently provided by receivers for 
acquiring seismic data in the time domain and transform means for transforming the 
data to the frequency-space domain. The transform means may be incorporated in the 
apparatus 1 1 or they may be additional to it. 

The apparatus 1 1 comprises a programmable data processor 12 with a program memory 
13, for instance in the form of a read only memory (ROM), storing a program for 
controlling the data processor 12 to process seismic data by a method of the invention. 
The apparatus further comprises non-volatile read/write memory 14 for storing, for 
example, any data which must be retained in the absence of a power supply. A 
"working" or "scratch pad" memory for the data processor is provided by a random 
access memory RAM 15. An input device 16 is provided, for instance for receiving 
user commands and data. One or more output devices 17 are provided, for instance, for 
displaying information relating to the progress and result of the processing. The output 
device(s) may be, for example, a printer, a visual display unit, or an output memory. 

Sets of seismic data for processing may be supplied via the input device 16 or may 
optionally be provided by a machine-readable data store 18. 

The results of the processing may be output via the output device 17 or may be stored. 

The program for operating the system and for performing the method described 
hereinbefore is stored in the program memory 13, which may be embodied as a 
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semiconductor memory, for instance of the well known ROM type. However, the 
program may well be stored in any other suitable storage medium, such as a magnetic 
data carrier 13a (such as a "floppy disk") or a CD-ROM 13b. 
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CLAIMS: 

1 . A method of processing seismic data, the data being in the frequency-space 
domain, the method comprising: 

determining a prediction filter from the seismic data at a first frequency at which 
swell noise is not present; and 

applying the prediction filter to seismic data at a second frequency at which 
swell noise is present, thereby to attenuate swell noise in the seismic data at the second 
frequency. 

2. A method as claimed in claim 1 and further comprising the steps of: 
determining a prediction filter from the seismic data at a third frequency at 

which swell noise is not present, the third frequency not being equal to the first 
frequency; and 

applying the prediction filter to seismic data at a fourth frequency at which swell 
noise is present, the fourth frequency not being equal to the second frequency, thereby 
to attenuate swell noise in the seismic data at the fourth frequency. 

3. A method as claimed in claim 2 wherein the difference between the first 
frequency and the second frequency is equal to the difference between the third 
frequency and the fourth frequency. 

4. A method as claimed in claim 1 , 2 or 3 and comprising, for each frequency in 
the range for which swell noise is present in the seismic data, determining an associated 
prediction filter from the seismic data at a frequency at which swell noise is not present. 

5. A method as claimed in claim 4 wherein the frequency range for which swell 
noise is present in the seismic data extends from a frequency to a frequency f H , and 
the associated prediction filter for frequency is determined from seismic data at a 

frequency fi-L+H+i- 
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6. A method of marine seismic surveying comprising actuating a source to emit 
acoustic energy; acquiring seismic data in the frequency-space domain; and processing 
the data according to a method as defined in any preceding claim. 

7. A method as claimed in claim 6 wherein the step of acquiring seismic data in the 

frequency-space domain comprises acquiring seismic data in the. time domain and 

! 

transforming the data to the frequency-space domain. 

8. A method as claimed in any preceding claim and comprising the further step of 
transforming the processed data to the time domain. 

9. An apparatus for processing seismic data, the data being in the frequency-space 
domain, the apparatus comprising: 

means for determining a prediction filter from the seismic data at a first 
frequency at which swell noise is not present; and 

means for applying the prediction filter to seismic data at a second frequency at 
which swell noise is present, thereby to attenuate swell noise in the seismic data at the 
second frequency. 

10. An apparatus as claimed in claim 9 and further comprising: 

means for determining a prediction filter from the seismic data at a third 
frequency at which swell noise is not present, the third frequency not being equal to the 
first frequency; and 

means for applying the prediction filter to seismic data at a fourth frequency at 
which swell noise is present, the fourth frequency not being equal to the second 
frequency, thereby to attenuate swell noise in the seismic data at the fourth frequency. 

11. An apparatus as claimed in claim 9 or 10 and adapted to determine, for each 
frequency in the range for which swell noise is present in the seismic data, an associated 
prediction filter from the seismic data at a frequency at which swell noise is not present. 
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12. An apparatus as claimed in claim 11 wherein the frequency range for which 
swell noise is present in the seismic data extends from a frequency f L to a frequency f H , 
and the apparatus determines, in use, an associated prediction filter for frequency fj from 
seismic data at a frequency fi-L+H+i- 

13. An apparatus as claimed in any of claims 9 to 12 and comprising a 
programmable data processor. 

14. A seismic surveying arrangement comprising: a source of acoustic energy; 
means for acquiring seismic data in the frequency-space domain; and an apparatus as 
defined in any of claims 9 to 13. 

15. A seismic surveying arrangement as claimed in claim 14 wherein the means for 
acquiring seismic data in the frequency-space domain comprises one or more receivers 
for acquiring seismic data in the time domain and transform means for transforming the 
data to the frequency-space domain. 

16. A storage medium comprising a program for a data processor of an apparatus as 
defined in claim 13. 

17. A storage medium containing a program for controlling a data processor to 
perform a method as defined in any of claim 1 to 5. 
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